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Abstract 

We report on the measurement of various flow harmonics, v„, with multi-particle cumulants, and 
present the results from a study of the inter-correlation among different order symmetry planes 
via multi-particle mixed harmonic correlations. This provides comprehensive experimental 
information on the fluctuating event-by-event shape of the initial conditions, which is currently 
among the main sources of large theoretical uncertainties in describing the evolution of the sys- 
tem created in heavy-ion collisions. 



1. Introduction 

The properties of the produced matter in heavy-ion collisions can be experimentally stud- 
ied by measuring the azimuthal anisotropy in the momentum distribution of the produced par- 
ticles UlEl. Quantified by the anisotropic flow coefficients, v„, and corresponding symmetry 
planes, *P„, such anisotropy is expected to reflect the shape and fluctuations in the initial energy 
density of the collision. Fluctuations in initial geometry contribute to the even harmonic flow 
and give rise to non-zero odd harmonic flow [3|. Anisotropic flow harmonics are calculated as 
an average v„ = (cos[h(i^- v F„)]), where <p's are azimuthal angles of all reconstructed particles in 
an event. 

The contribution from fluctuations implies that ^v*^ is not the same as (v„) k , and the sym- 
metry planes T*,, are in general different. In these proceedings we will search for answers to 
the following two questions: a) What is the underlying probability density function (p.d.f.) of 
flow fluctuations which uniquely determine all moments ^v^?; b) What is the relation between 
different symmetry planes ^P,,? 

The cornerstone of our approach is the following relation for multi-particle correlations: 

<cos(«i^i + ■ • ■ + n k ip k )) = v„, ■ • ■ v„ k cos(n l x l J „ 1 + ■ • • + n^^) , (1) 

derived first in |4|. The brackets on the LHS indicate the single-event average over all distinct 
tuples of k particles. The last term in Eq. ([T]i on the RHS quantifies the inter-correlations among 
different order symmetry planes ^P,,. 
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One may wonder to what extent the underlying p.d.f. of v„ fluctuations can be uniquely 
reconstructed from the measured moments (vty, or better from the measured multi-particle cu- 
mulants (which are less sensitive to unwanted few particle non-flow correlations). Whether or 
not the equivalence "p.d.f. o moments o cumulants" holds true was subject to debate for a long 
time among mathematicians. It turns out that the 2 nd equivalence is always true, whilst the 1 st one 
is not true in general. However, given the functional form of p.d.f., one can easily test whether it 
is uniquely determined via its moments by checking Krein-Lin conditions [5 1. The most popular 
candidate on the theoretical frontline for the p.d.f. is a Bessel-Gaussian function, which has the 
remarkable property that all its multi-particle cumulants are the same |6|. The Bessel-Gaussian 
function satisfies the Krein-Lin conditions, which means that it is a unique p.d.f. for which all 
multi-particle cumulants are the same. 



2. Analysis details 

We have used a multi-particle azimuthal correlation technique Q. The analysis was per- 
formed over the dataset comprising about 30 million Pb-Pb collisions at y/s^N = 2.76 TeV both 
from the 2010 and 201 1 data taking periods. Data were taken with the ALICE detector (8). Se- 
lected tracks of charged particles were reconstructed using the ALICE Time Projection Chamber 
(TPC) |9|, while for centrality determination we have used the VZERO detectors. Only events 
with reconstructed primary vertex within 10 cm along the beam line from the nominal interaction 
point were selected. Charged particle tracks were required to have at least 70 reconstructed space 
points in the TPC and a (x 2 ^J per TPC cluster < 4. The maximal distance of closest approach of 
reconstructed tracks to primary vertex, both in longitudinal and transverse direction, was set to 
3.0 cm. Event and track selection parameters were varied to estimate the systematic error of the 
measurement. 
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Figure 1: (color online) Centrality dependence of vi, V2 and 
V3 estimated with multi-particle cumulants. 



Figure 2: (color online) Transverse momentum dependence 
of v'2 and V3 estimated with 4-particle cumulants. 



3. Results 

Figure[T]shows results for vi, V2 and V3 estimated with multi-particle cumulants (such results 
for V2 were already available in the analysis for iflOll ). We observe that for all centralities and 
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harmonics v„{4) v„{6). This implies that the underlying p.d.f. of flow fluctuations must have 
non-negligible 3 rd /higher moments when compared to the l st /2 nd moment. A Bessel-Gaussian 
p.d.f. is consistent with these experimental findings. 

Figure [2] shows V2{4) and V3{4) measured as a function of transverse momentum for three 
different centralities, in the range 0.2 < pj < 5 GeV/c. The coefficient V3{4) increases steadily 
up to pi ~ 3 GeV/c and then saturates. Similar behavior was found for V2{4), which suggests 
a similar origin of both type of anisotropics. Compared to v'2{4), v^,{4}(pj) changes little with 
centrality. This is consistent with the widely accepted picture in which V3 originates mainly from 
anisotropics stemming from event-by-event fluctuations of the initial positions of participating 
nucleons, while v>2 has an additional large contribution from almond shape geometry of the nuclei 
overlap region in non-central collisions. 
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, <cos(3<p +3<f> z -2ip 3 -2<p 4 -2c|> s )> (PRL 107 (2011) 032301) 
<cos(2tp ] +2(p z -2(p 3 -(p 4 -(p 5 ))-2(cos(2(p ] -2tp 2 ))(cos(2(() i -(p 2 -tp 3 )) 
; [^] (cos(3(p ] +2^-2ip 3 -2(p j -ip s ))-2{cos(2(p ] -2(p^))(cos(3(p i -2(P 2 -{p 3 )) 
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Figure 3: (color online) Transverse momentum dependence Figure 4: (color online) Centrality dependence of 5-particle 
of 3-particle correlator. cumulants. 

Figure [3]presents the 3-particle correlator (cos(^„— 3^j+2 1 F2)), which was proposed first by 
Teaney and Yan in ifTTIl . By making comparison with Fig. 15 from [ 1 1], we conclude that ideal 
hydrodynamic model calculations with Glauber initial conditions can only qualitatively describe 
the data at low p T < 2 GeV/c, while the agreement is rather poor at higher p T . 

Figure|4]presents our results obtained with 5-particle cumulants, which are sensitive to inter- 
correlation among different order symmetry planes *?/„, in particular to VjV, cos[6( v F3- v f , 2)] (green 
markers), — v|vj cos[2( v F2- v i'i)] (red markers) and -V3v|vi cos [3*1/3- 2*^2 -^i] (blue markers). 
Results on 3- and 5-particle correlations show the existence of strong genuine three plane cor- 
relation (cosO?! -3 ( f , 3+2 v r'2)) among symmetry planes Ti, ^2 and ¥3 in mid-central collisions, 
although the two plane correlation (cos[6(*f2 - ¥3)]) between symmetry planes ¥2 and ¥3 is 
negligible in the same centralitites. 



4. Summary 

We have reported in this contribution detailed results of anisotropic flow analysis which can 
further improve our understanding of the properties of the system produced in relativistic heavy- 
ion collisions. The probability density function of event-by-event flow fluctuations was studied 
with multi-particle cumulants, and was found to be consistent with a Bessel-Gaussian distri- 
bution. The transverse momentum dependence of triangular flow measured with 4-particle cu- 
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mulants was compared to analogous previous measurements for elliptic flow, and it was demon- 
strated that these two observables reflect different physics, and can be used independently to con- 
strain outputs of theoretical models. Multi-particle cumulants in mixed harmonics were reported 
as well, and it was shown that the theoretical framework [ 1 1 1 which relies on ideal hydrodynamic 
model calculations with the Glauber initial conditions can neither quantitatively reproduce the 
measurements at low pj < 2 GeV/c, nor even qualitatively at higher pj. 
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